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We present a comprehensive mean-field calculation of the Schiff moment of the nucleus ^^®Ra, 
the quantity which determines the static electric dipole moment of the corresponding atom if time- 
reversal (T) invariance is violated in the nucleus. The calculation breaks all possible intrinsic 
symmetries of the nuclear mean field and includes, in particular, both exchange and direct terms 
from the full finite-range T-violating nucleon-nucleon interaction, and the effects of short-range cor¬ 
relations. The resulting Schiff moment, which depends on three unknown T-violating pion-nucleon 
coupling constants, is much larger than in ^®®Hg, the isotope with the best current experimental 
limit on its atomic electric-dipole moment. 


The Standard Model of particle physics violates time- 
reversal (T) invariance, but apparently only through a 
single phase in the Cabibo-Kobayashi-Maskawa matrix 
that mixes quark flavors. The resulting T violation in 
flavor-conserving observables is therefore very weak and 
static electric-dipole moments (EDMs) of neutrons, elec¬ 
trons, or atoms, all of which are nonzero if T is violated, 
have never been observed. Standard-Model T violation 
is also too weak to account for the baryon asymmetry 
of the universe, which must come from as yet undiscov¬ 
ered physics. Happily, most theories of what lies be¬ 
yond the Standard Model contain enough phases that 
flavor-conserving T violation will be unsuppressed. Cur¬ 
rent levels of sensitivity in EDM experiments are already 
sufficient to rule out or pressure many extra-Standard 
models, and it seems quite possible that with slightly 
improved sensitivity, new T-violating physics will be dis¬ 
covered. 

Some of the tightest constraints on T violation come 
from atomic EDM experiments. The best of these at 
present is an experiment [1] with ^®®Hg, but it has 
become clear recently [2-4] that atoms with octupole- 
deformed nuclei are potentially more sensitive than Hg. 
The primary reason is that given any T violation in 
the nucleon-nucleon interaction, an asymmetric nuclear 
shape and an associated parity doubling create a collec¬ 
tive “Schiff” moment, a kind of radially weighted dipole 
moment (see below). [Earlier papers on T violation in 
polar diatomic molecules [5] and parity violation in fis¬ 
sion [6] were the first to consider the effects of asymmetric 
shapes on fundamental symmetries. Reference [7] drew 
early attention to the importance for EDMs of low-lying 
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nuclear states with the same angular momentum as the 
ground state but opposite parity.] Because of screening 
by atomic electrons, the Schiff moment, rather than the 
nuclear EDM, is the quantity that directly induces an 
atomic EDM (at least in lowest order; see Ref. [8]). In 
nuclei with symmetric shapes, a collective contribution to 
the Schiff moment develops only in fluctuations around 
that shape [4, 9]. 

In this paper we calculate the Schiff moment of ^^®Ra, 
or more precisely its dependence on any T-violating ttNN 
couplings, in a mean-field theory that allows us to break 
all possible symmetries, consider a variety of phenomeno¬ 
logically successful strong (Skyrme) interactions, implic¬ 
itly include the RPA polarization of the even-even core by 
the valence neutron, treat both the direct and exchange 
parts of the full pion-mediated interaction responsible for 
creating the Schiff moment, and include short-range two- 
body correlations between nucleons that modify the ef¬ 
fects of this T-violating interaction. Though further re¬ 
finements are possible, they will probably have to include 
correlations beyond mean-field theory and/or careful and 
systematic work on Skyrme functionals; the results pre¬ 
sented here will not be easy to supersede. 

Simpler calculations of Schiff moments have been at¬ 
tempted before. Ref. [10] applied an independent-particle 
model in ^®®Hg, ^^®Xe, and other symmetrically de¬ 
formed or spherical isotopes. Refs. [11, 12] carried out 
a much more sophisticated RPA-based calculation in 
i99Hg; it’g main drawback was the use of a single phe¬ 
nomenological interaction that made it difficult to esti¬ 
mate uncertainty. Ref. [3] made estimates in a particle- 
rotor model of the enhancement due to octupole deforma¬ 
tion, and in Ref. [13] we applied a preliminary version of 
our technique to ^^®Ra, an experiment on which is in the 
works [14]. That paper, however, assumed the range of 
the T-violating interaction to be zero, an especially bad 
approximation for exchange matrix elements, and was 
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unable (obviously) to examine the effects of short-range 
NN correlations. 

We briefly review some definitions and ideas. The 
Schiff moment is given by 

5 . = E + C.C., 
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where j'l'o) is the member of the ground-state multiplet 
with Jj = J = 1/2 (positive parity), the sum is over 


excited states, and Sz is the operator 

= ( 2 ) 
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with the sum here over protons, and the mean-square 
charge radius. The operator Vpx in Eq. (1) is the T- (and 
parity-) violating nucleon-nucleon interaction mediated 
by the pion [7, 15]: 


Vpriri -r2) = - 


gm; 


{(^ 


I \Cri - cr2) ■ (ri - r 2 ) gon ■ T2- ^ {tiz + T2z) + g2{iTizT2z - Ti ■ T 2 ) 

oTrmjv I L i 


( 3 ) 


-y (o-i -f (T2) ■ (ri - r2) (tiz - T2z)| 


'I exp{-m^\ri - r2\) 


m„\ri - r2 


1-f 


m-^\ri - r2 


where arrows denote isovector operators, Tz is -|-1 for neu¬ 
trons, ruAT is the nucleon mass, and (in this equation 
only) we use the convention h = c = 1. The g's are the 
unknown isoscalar, isovector, and isotensor T-violating 
pion-nucleon coupling constants, and g is the usual strong 
ttNN coupling constant. 

The asymmetric shape of ^^®Ra implies parity dou¬ 
bling (see e.g. Ref. [16]), i.e. the existence of a very low- 
energy |l/2“) state, in this case 55keV [17] above the 
ground state Id^o) = 11/2+), that dominates the sum in 
Eq. (1) because of the corresponding small denominator. 
With the approximation that the shape deformation is 
rigid, the ground state and its negative-parity partner in 
octupole-deformed nucleus are projections onto good par¬ 
ity and angular momentum of the same “intrinsic state” 
(see Fig. 1), which represents the wave function of the nu¬ 
cleus in its own body-fixed frame with the total angular 
momentum aligned along the symmetry axis. Equation 
(1) then reduces to [3] 
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where the brackets indicate expectation values in the in¬ 
trinsic state. Using Eq. (3) for Vpt, we can express the 
dependence of the Schiff moment on the undetermined 
T-violating ttNN vertices as 


S = aoggo + aiggi + a2gg2 ■ (5) 


The coefficients Oi, which are the result of the calculation, 
have units e fm^. 

The octupole deformation enhances (Sz), the first fac¬ 
tor in Eq. (4), making it collective, robust, and straight¬ 
forward to calculate with an error of a factor of two or 
less. The interaction expectation value {Vpt) is harder 
to estimate because it is sensitive to the nuclear spin 



FIG. 1: (color online). Shape of the microscopically calcu¬ 
lated [13] mass distribution in ^^®Ra, represented here by the 
surface of a uniform body that has the same multipole mo¬ 
ments Qao for A=0... 4 as our calculated density. 


distribution, which depends on delicate correlations near 
the Fermi surface. Our calculation allows the breaking 
of Kramers degeneracy in the intrinsic frame and, conse¬ 
quently, spin polarization. 

To evaluate (Vpt) we constructed a new version of the 
code HFODD (v2.14e) [18, 19]. The code uses a triax- 
ial harmonic-oscillator basis and Gaussian integration to 
solve self-consistent mean-field equations for zero-range 
Skyrme interactions. Evaluating matrix elements of the 
finite-range interaction (3) is much harder numerically, 
but efficient techniques have already been developed [20] 
for Gaussian interactions, which are separable in three 
Gartesian directions. The spatial dependence in Eq. (3) is 
different, the derivative of a Yukawa function, and we also 
include short-range correlations between nucleons (which 
the mean-field does not capture) by multiplying the in- 
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teraction by the square of a correlation function [21] that 
cuts off the two-nucleon wave functions below a relative 
distance of about a fermi: 

/(^) = l_e-ii’-=(l-.68r2) , (6) 


with r = jri — r 2 | in fermis and the coefficients of in 
fm“^. The resulting product looks very different from a 
Gaussian, but we were able to reproduce it quite accu¬ 
rately (see Fig. 2) with the sum of four Gaussians: 


J 
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where Ott = 0.7045 fm~^ is the pion mass in inverse fermis 
and the numbers in the fit all have units fm“^. The extra 
factor of r not included in Eq. (7) (i.e. the factor ri — r 2 
in Eq. (3)) is treated separately. 

HFODD works with any Skyrme energy functional. In 
the context of the present study the best is SkO' [22, 23]. 
The “time-even” terms in this interaction, which act in 
even-nucleus ground states, were fit with special atten¬ 
tion to nuclei around ^°®Pb and to spin-orbit splitting. 
The “time-odd” terms responsible for core-polarization 
in an odd nucleus were adjusted in Ref. [23] to reproduce 
Gamow-Teller resonances, resulting in an effective Lan¬ 
dau parameter = 1.2. (The isoscalar parameter was 
set to go = 0.4, following common practice.) For com¬ 
parison, we also carry out the calculation with the older 
parameterizations Sill, SkM*, and SLy4 (with time-odd 
terms that are not fixed by gauge invariance neglected, 
and then again, with the simplest time-odd terms mod¬ 
ified so that the Landau parameters have the same val¬ 
ues as in SkO') but have the most confidence in SkO'. 
Ref. [13] presented predictions by these functionals for 



FIG. 2: (color online). The function g{r) in Eq. (7) multiplied 
by r® (solid line), the Gaussian fit multiplied by r® (dashed 
line), and g{r)/f{r)^, the radial T-odd interaction without 
short-range correlations (dot-dashed line). The factor r® is to 
account for the volume element and the additional factor of 
r = ri — r2 in Eq. (3). 


-f 0.53 0.11 0.004 (7) 


the binding energies, separation energies, intrinsic dipole 
moments, and spin-orbit splittings in the even Ra iso¬ 
topes. SkO' and Sill seemed to do the best job. 

TABLE I: Coefficients of ggi, in units of efm®, in the expres¬ 
sion Eq. (5) for the Schiff moment of ^^®Ra, calculated with 
the SkO' Skyrme interaction. The abbreviation “src” stands 
for “short-range correlations”. 



ao 

ai 

a2 

zero-range (direct only) 

-5.1 

10.4 

-10.1 

finite-range (direct only) 

-1.9 

6.3 

-3.8 

finite-range -|- src (direct only) 

-1.7 

6.0 

-3.5 

finite-range -|- src (direct -|- exchange) 

-1.5 

6.0 

-4.0 


Table I shows the calculated values, with SkO', of the 
three coefficients Ui at several levels of approximation. 
The finite range reduces the direct matrix elements of 
the interaction (and the corresponding a^) from the zero- 
range limit [13] significantly. The exchange terms are re¬ 
duced much more, so that they are always smaller than 
the direct terms. The effects of the short-range corre¬ 
lations, which also reduce the coefficients, are relatively 
small as well but non-negligible. Finally, as pointed out 
in Ref. [8], relativistic effects in electron wave functions 
correct the effects of the Schiff moment; the authors sum¬ 
marize the corrections in a quantity they call the “local 
nuclear dipole moment”. Our local dipole moment in 
^^®Ra is 82% of the Schiff moment for the Si /2 — P 1/2 
atomic transition and 87% for the Si /2 —P 3/2 transition. 


TABLE II: Same as the last line in Table I but for the Sill, 
SkM*, and SLy4 interactions. 



ao 

ai 

a2 

Sill 

-1.0 

7.0 

-3.9 

SkM* 

-4.7 

21.5 

-11.0 

SLy4 

-3.0 

16.9 

1 

00 

bo 


The coefficients produced by the forces we favor less, 
with the time-odd terms adjusted as mentioned above, 
are presented in Tab. II; Sill produces numbers similar 
to those of SkO', while the other two forces give numbers 
that are larger by factors of two or three, whether or not 
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the Landau parameters are adjusted (i.e., the effects of 
adjusting those parameters seem to be fairly small). 

What is the uncertainty in our numbers? The mean- 
field omits correlations that could have some effect on the 
result; those could be explored by refining the calculation 
through angular-momentum and parity projection, i.e., 
the restoration of symmetries broken by the mean field. 
In addition, an optimal Skyrme functional has yet to be 
identified. Those we tested give results that differ from 
the SkO' numbers by factors two or three, as mentioned 
above. But some low-order terms in the T-odd part of the 
Skyrme functional are never used even in SkO', because 
they have never been fit. Ref. [23] constrained some com¬ 
binations of those terms but others were set to zero for 
lack of sufficient Gamow-Teller data in spherical nuclei. 
One might imagine trying to fit in deformed nuclei, or 
looking at spin-strength distributions with different total 
angular momentum and parity; the 0“ channel would be 
particularly useful because those are the quantum num¬ 
bers of VpT- At the same time, it would probably help 
to explicitly study the sensitivity of the Schiff moments 
to changes in the various Skyrme parameters, both in 
the time-odd and time-even sectors. With enough work 
on all these fronts we could give a firmer estimate of the 
uncertainty than our current guess: a factor of two or 
three. 

If we accept our current results as reasonably accurate, 
we are in a position to quantify the advantages of ^^®Ra 
for an EDM measurement. A recent RPA calculation 
[II, 12] of the Ui for ^^^Hg gives 

i®®Hg : ao = 0.0004 , ui = 0.055 , 02 = 0.009. (8) 


Our numbers are more than two orders of magnitude 
larger, particularly in the isoscalar channel (the relevant 
channel if T-violation is caused by a nonzero QCD pa¬ 
rameter 0), where the Hg number is anomalously small. 
Atomic physics enhances any EDM in Ra by another fac¬ 
tor of 3 over that in Hg [24], so if the Ra EDM can 
be measured even 1/lOOth as accurately as that of Hg, 
the sensitivity to nuclear T violation will be significantly 
greater. 

In conclusion, we have evaluated the Schiff moment 
of ^^®Ra in a completely symmetry-breaking mean-field 
approach, including in particular the finite-range matrix 
elements of the T-violating nucleon-nucleon interaction. 
The results indicate that EDM experiments in this sys¬ 
tem are very promising. The remaining uncertainties of a 
factor of two or three are related primarily to deficiencies 
in nuclear effective interactions, which can be removed 
but not easily. 
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